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We study the angular process related to random walks in the Euclidean and in the non-Euclidean 
^~ 5 space where steps are Cauchy distributed. 

|y\ This leads to different types of non-linear transformations of Cauchy random variables which 

£N| preserve the Cauchy density. We give the explicit form of these distributions for all combinations of 

the scale and the location parameters. 

Continued fractions involving Cauchy random variables are analyzed. It is shown that the n-stage 
random variables are still Cauchy distributed with parameters related to Fibonacci numbers. This 
permits us to show the convergence in distribution of the sequence to the golden ratio. 

Keywords: hyperbolic trigonometry, arcsine law, continued fractions, Fibonacci numbers, non-linear trans- 
formations of random variables. 

> 

^h 1 Introduction 

ON 

We consider a particle starting from the origin O of K 2 which takes initially a horizontal step of length 1 

r- — , and a vertical one, say Ci, with a standard Cauchy distribution. It reaches therefore the position (1, Ci). 

The line l\ joining the origin with (l,Ci) forms a random angle Oi with the horizontal axis (See Figure 

On li the traveller repeats the same movement with a step of unit length (either forward or backward) 
£> along li and a standard Cauchy distributed step, say C2, on the line orthogonal to l\. The right triangle 

obtained with the last two displacements has an hypothenuse belonging to the line li with random 
inclination O2 on l\. 

After n steps the sequence of random angles ©i, • • • , ©„ describes the rotation of the moving particle 
around the starting point, their partial sums describe an angular random walk which can be written as 

n 

S„ = @i-\ h 9„ = ^2 arctanCj (1.1) 
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Figure 1: The angular process in the Euclidean plane. By Cr ai b) we indicate the j-th random displacement 
with Cauchy distribution possessing scale parameter dj and location parameter bj. 



where Cj are independent standard Cauchy random variables. If the random steps of the planar random 
walk above were independent Cauchy random variables with scale parameter a,j and location parameter 
bj then the process ( 1.1 ) must be a little bit modified and rewritten as 



S n = 6i 



Qn 



y arctan Cj , 

3=1 



(1.2) 



where C, ~ C, 



( a j,bj)- 



The model ( 1.2 1 can be extended also to the case where the first step has length dj 



and the second one is Cauchy distributed with scale parameter Qj and position parameter bj (see Figure 



0. 



then 



tan 0j = C i 



The same random walk can be generated if the two orthogonal steps, at each displacement, are 
represented by two independent Gaussian random variables Xj and Yj. In this case, for each right 
triangle, we can write 

Xi 

tan 6,- = T7-, 

Ya 



If Xj and Yj are two standard independent Gaussian random variables then tan Qj = y 2 - possesses stan- 
dard Cauchy distribution and we get the model in ( 1.1 ). The model ( 1.2 ) can be obtained by considering 



orthogonal Gaussian steps with different variances and in this case the scale parameter of the random 



variables Cj is the ratio a,j 



1° 



The model (1.1) describing the angular random process has an hyperbolic counterpart. We consider 
a particle starting from the origin O of the Poincare half-plane H^ = {(x,y) : y > 0}. At the j-th 
displacement, j — 1,2..., the particle makes two steps of random hyperbolic length rjj and fjj on two 
orthogonal geodesic lines. The j-th displacement leads to a right triangle Tj with sides of length r\j and 
r)j and random acute angles Qj and Qj. In each triangle Tj the first step is taken on the extension of 
the hypotenuse of the triangle Tj-\ (see Figure |2J|. From hyperbolic trigonometry (for basic results on 
hyperbolic geometry see, for example, Faber [S]) we have that 



sin0j = 



sinh 



>li 



cosh rjj cosh fjj — 1 



cos Qj = 



sinh rjj cosh 



>li 



cosh rjj cosh fjj — 1 



From the above expressions we have that 



tan©i 



t&iahfjj 
sinh rjj 



If we take independent random hyperbolic displacements rjj and m such that the random variables 
E-i -■- jf " u zr are standard Cauchy distributed then Qj = arctan Cj. If the triangles Tj were isosceles 
and the angle Qj 6 [—?»?] so that in this case the Cauchy distribution cannot be 



3 sinh rij 

then tanO, = — \ — 

J cosh r\j 

attributed to tan@,-. 





Figure 2: The angular random process in the Poincare half-plane. 

In the model described here the random steps (and therefore the random angular windings Qj) are 
independent. If we consider the model of papers [2] and [3J, where the displacements are taken orthogo- 
nally to the geodesic lines joining the origin O of H^_ with the positions occupied at deviation instants, 
the angular displacements Qj must be such that 



sin Qj 



sinh 



Vj 



/ 3 

= sinh rjj cos arctan I I cosh r\ r 

1 + llr=l cosh2 Vr V 



r=\ 



and therefore involve dependent random variables. 

For the area of the random hyperbolic triangle Tj we note that 



area(7}) = 7 L-Q j -Q j = 7 L- 



arctan 



TV 

2 



/ COth 77, 

arctan . 

\ sinh rjj 



I tanh rjj 
\ sinh rjj 

coth fjj 

sinh rjj 



arctan 



/ tanh rjj 



arcotan 



\ sinh fjj 
( coth rjj coth fjj 
\ sinh fjj sinh rjj 



Since each acute angle inside Tj is linked to both sides of the triangle, the analysis of the random process 
A n = Ylj—! area(Tj) is much more complicated and we drop it. 



Let C, - C, 



(a, ,bj 



7 = 1,2... be independent Cauchy random variables where a,j is the scale parameter 



and bj is the location parameter. In the study of the angular random walk ( 1.1 1 and ( 1.2 ) we must analyze 
the distribution of the following non-linear transformations of Cauchy random variables: 



U 



C1+C2 
1 - C!C 2 



(1.3) 



since 



arctan C\ + arctan Ci = arctan 



Ci+C 2 



1 - CiC 2 

Since the Sixties a wide class of non-linear transformations of Cauchy random variables has been 
considered. Williams [TU], Knight [5] and Letac [7] proved that transformations of the form 



ef(x) = kx + a — \^ 



Hi 



where e = ±1, a, 7$ € R and fc,/^ > preserve the Cauchy distribution. 

In particular, in Williams [TU] it is proved the following characterization for Cauchy random variables. 
The random variable X is a standard Cauchy if and only if (1 + bX)/(b — X) is a standard Cauchy for 
some constant b which is not the tangent of a rational multiple of n. 

Knight [6] asserts that a random variable X is of Cauchy type if and only if the random variable 
(aX + b)/(cX + d) is still of Cauchy type, whenever ad — be 7^ 0. 

Our problem is more strictly related to the results obtained by Pitman and Williams [8 . They proved 
that the standard Cauchy distribution is preserved under certain types of transformations represented 
by meromorphic functions whose poles are all real and simple. As a corollary they obtained that, if P 
and Q are two independent random variables uniformly distributed in (0, 7r), then the random variables 
X = tanP and Y = tanQ are standard Cauchy and 



We will show that the random variable (1.3 1 is endowed with Cauchy distribution in a much more 
general situation, namely when the random variables Cj, j = 1,2, have non-zero location parameters bj 
and scale parameters a,j. The scale and location parameters of U depend on both parameters a, and bj 
suitably combined. 

In particular, if bi = b- 2 = and a\ — a 2 = 1, then U is still distributed as a standard Cauchy 



distribution and therefore in (1.1) we have that 



S n = arctan ( 



Since also i is a standard Cauchy (for basic properties of Cauchy random variables see, for example, 
Chaumont and Yor [5] page 105), from ( |1.3| , a number of related random variables preserving the form 
of the Cauchy distribution can be considered. For example, the following random variables 



_ dC 2 + 1 _ 1-' 

£\ — — — — , /12 — -=; — 



also possess standard Cauchy distribution. 



If b\ = 62 = and the scale parameters ai,a2 are different, then (1.3 1 still preserves the Cauchy 
distribution but with scale parameter equal to 1 °?;^" a a 2 and location parameter equal to zero. This can be 
grasped by means of the following relationship 

arctan Ci + arctan C2 == arctan < C > , (1-4) 

11 + axa 2 I 



where C, ~ Ct a . ) ■ Result ( 1 .4 



By iterating the process (|1.4[) we arrive at the formula 



is illustrated in Figure 



V- _^ n. '£? __ / ^U a > + a i a ^ 3 



arctan C,- = arctan 



^ ' 3 1 i + Ei#i<»iOi 



which gives an insight into further extensions of the process outlined above. 

Much more complicated are the cases where the location parameters of the Cauchy distributions are 
different from zero. For the special case where C\ and C 2 are independent Cauchy such that C\ ~ C(i,6) 
and C2 ~ C(i,6)) the random variable (1.3 1 still possesses Cauchy density with scale parameter 2 b b 4 T^ 4 and 
location parameter ,4°. . 
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Figure 3: The figure shows that shooting a ray with inclination Or, uniformly distributed, against the line at 
distance a\ and then shooting a ray with a uniformly distributed angle O2 on the line at distance 0,2 is 
equivalent to shooting on the barrier at the distance a = /V^ 2 with a uniformly distributed angle O. 



and Co 



We have obtained the general distribution of (1.3 1 where C\ and Ci are independent Cauchy such 
that C\ ~ C(oi,6i 



C(a 2 ,b 2 ) an d also the distribution of 
u= jC 1 + 6C 2 



a — jiC\Ci 

for arbitrary positive real numbers /3, 7, (5 and a + f3 =/= 0. In particular, if Ci and Ci are independent 
standard Cauchy random variables then U is Cauchy with scale parameter equal to 
parameter equal to zero. 



-y+S 



and location 



In the last section we have examined continued fractions involving Cauchy random variables. In 
particular we have studied 

V„ = — x (1.5) 



and 



U n 



1 



1 



(1.6) 



which generalize the random variables V\ = y^ and U\ = . ,^2 ■ Continued fractions involving random 
variables have been analyzed in Chamayou and Letac [1] and more recently in Asci, Letac and Piccioni pQ. 
The random variable U\ has the arcsine distribution in [0, 1], while XJ t — tUi, with t > 0, has distribution 



Pr{C7 t E ds} = 



d.s 



TTy/s{t~ S) 



< s < t. 



For each n > 1, the random variables V n , are Cauchy distributed with scale parameter a n and position 
parameter b n that can be expressed in terms of Fibonacci numbers. This permits us to prove the mono- 
tonicity of a n and b n and that lim„_ i . 00 a n = and lim„_>. 00 b n — tfi — 1 where <fi = ^ I s tl ie golden 
ratio. Finally we obtain that the sequence of random variables 1 + V n and 1 + U n , n > 1, converges 
in distribution to the number <j) = 1+ ^ ■ This should be expected since it has the infinite fractional 
expansion 

l + \/5 . 1 



= 1 + 



1 



(1.7) 



which is related to (1.5) and (1.6 1 



2 Centered Cauchy random variables 

In this section we study the distribution of the following random variable 

7C1 + 5C 2 



U 



'l*-2 



(2.1) 



where Ci and C2 are independent, standard Cauchy. We assume, without restrictions, that /3, 7, 5 are 
non-negative real numbers all different from zero (because of the symmetry of Cj, j — 1, 2) and a + (3 ^ 0. 
In the next theorem we prove that U is still Cauchy distributed. 



7+5 
Q+/3 



Theorem 2.1. The random variable U in (2.1) possesses Cauchy distribution with scale parameter equal 
to 



and position parameter equal to zero. We can also restate the result in symbols as 

law 7 + ^ 



U 



a + /3 



Proof 1 
The density of the random variable U can be obtained by means of the transformation 



-fx+Sy 
a—j3xy ' 



with (x, y) £ M 2 and Jacobian equal to 



\J\ 



The joint density g = g(u, v) reads 



g(u,v) 



1 



aj + 8j3v 2 
(7 + j3uv) 2 

1 «7 + 5{3v 2 



TT 2 {1 + V 2 ) f au -Sv \ 2 {1 + PUV) 2 

+ I 7+0uv I 



(u,v) G 



1 



«7 + 5(3v 2 



7r 2 (l + v 2 ) v 2 (/3 2 u 2 + S 2 ) + 2uv(f3~f - da) + a 2 u 2 + 7 2 
and must be integrated with respect to v in order to obtain the distribution of U. Therefore 



[a- p<Li<C2 



du 


pOO 


7T 2 „ 


— 00 


du 


/•OO 


7T 2 „ 


— 00 



1 



aj + 8(3v 2 



1 + v 2 v 2 (u 2 /3 2 + 5 2 ) + 2uw(/?7 - 5a) + u 2 a 2 + 7 2 
Av + B Cv + D 



where 



r a 



l + v 2 v 2 {u 2 I3 2 + 5 2 ) + 2vu(/3j - 5a) + u 2 a 2 + 7 2 

2u(P*f-a8)(p8-a'y) 

[„2( Q _ /3 )2 + ( 7 _ (5) 2][ u 2( a + /3 )2 + ( 7+( 5)2], 



dv (2.2) 
dv 



C = -A(u 2 f3 2 + 5 2 ), 
D 



[« 2 (a-^) 2 + ( 7 -<5) 2 ][« 2 (a+/3) 2 + (7+<5) 2 ] 



We start by evaluating the first part of the integral (2.2 1 

Av 



Cv 



1 + v 2 v 2 (u 2 (3 2 + S 2 ) + 2vu(P"f — 5a) + u 2 a 2 + 7 2 



dv 



A < x 

2 

A 



— lim log 

2 d— >oo,c— y— oo 

+Au((3j ~ 5a) 
Au{ji~i - 5 a) 



2v 2v{u 2 P 2 + 5 2 ) ± 2m(^ 7 - 5a) 

l + v 2 v 2 (u 2 j3 2 + 5 2 ) + 2vu(/3j - 5a) + u 2 a 2 + 7 2 

1 + v 2 



dv 



2 d^oo, c ^-co b \v 2 (u 2 /3 2 + 5 2 ) + 2vu(f3j - 5a) + u 2 a 2 + 7 2 

1 
v 2 {u 2 j3 2 + 5 2 ) + 2vu(/3j - 5a) + u 2 a 2 + -f 2 
1 



dv 



_ 00 v 2 {u 2 f3 2 + <5 2 ) + 2vu(/3j — 5a) + u 2 a 2 + y 



:dv 



Au(/3j - 5a) 



u 2 a(3 + jd 



where the last integral is obtained by means of the change of variable 



/ „ „„ ~ uifll — 8a) / 

vJu 2 P 2 + 8 2 + ) ' -. = z\ ■ 

^u 2 p 2 +5 2 V 



= z\ I u 2 a 2 + 7 2 



u 2 (/3~/ — 5a) 2 u a/3 + y5 

u 2 f3 2 + 5 2 ~ Z y/uipz + S 2 ' 



In view of result (2.3 1 and inserting the values of A, B and D we have that 

'Av + B 



dv 

~~~2" 

' — oo 

du 1 



Cv + D 



7r u 2 a(3 + 7<5 
du 

IT 



1 + v 2 v 2 (u 2 (3 2 + 5 2 ) + 2vu(/3"f — 5a) + u 2 a 2 + j 2 
[Au{Py - aS) + B(u 2 a(5 + yd) + D] 



dv 



( 7 <5 + M 2 a/3)(/3 7 - a8)[u 2 {fi 2 - a 2 ) + ( 7 2 - <5 2 )] 



(w 2 a/3 + -/S)[u 2 (a - /3) 2 + ( 7 - 5) 2 }[u 2 (a + /3) 2 + ( 7 + £) 2 
( 7 <5 + M 2 afi)(a 7 - /3<5)[ M 2 (a 2 - /3 2 ) + ( 7 2 - 5 2 )} 
+ (u 2 a(3 + j5)[u 2 (a - (3) 2 + (7 - 5) 2 }[u 2 (a + /3) 2 + (7 + <5) 2 ] 
(to ^ 2 (/3 2 - a 2 ){P - a)( 7 + J) + ( 7 2 - ^ 2 )(a + /3)( 7 - 5) 
■k [ M 2 (a-/3) 2 + ( 7 -(5) 2 ][u 2 (a + /3) 2 + ( 7 + 5) 2 ] 
du (a + /3)( 7 + <5)[u 2 (a-/3) 2 + (7-^) 2 ] 
7T [u 2 (a-/3) 2 + (7-,5) 2 ][M 2 (a + /3) 2 + (7 + <5) 2 ] 
du (a + )3)(7 + (J) 
~~ 7T w 2 (a + /3) 2 + ( 7 + <5) 2 ' 

Another approach is based on the conditional characteristic function, 

1 f°° jsf T"+ a " 1 du 

7T J_ 00 1 + ?i^ 

I /°° c ^ n/ ^ ^ ^ a 1 + 5Pv 2 

IT 



E{e^ u \C 2 =v} 



-du> 



(2.3) 



v 2 (/3 2 w 2 + 5 2 ) + 2vw(f3"/ — 5a) + u> 2 a 2 + 7 ; 
where w = 2 «„." ■ The inverse Fourier transform gives the conditional density g(w\v) and thus we arrive 



ca—fHuv ' 



again at the integral (2.2 1 



Remark 2.1. A special case implied by Theorem |2. 1 1 concerns the random variable 



V 



1 - C X C 2 



where Cj ~ C(o,o)) with j = 1,2, are independent Cauchy random variables with location parameter 
equal to zero and scale parameter aj > 0. If we choose 7 = ai, 5 — 0,2, a — 1 and /9 = 0102, we conclude 



that 

y law aiCi+02Ci law 0,\ + a 2 _ 

l-aia 2 CiC 2 1 + aia 2 

In the same way we obtain that 

c- _ 7C1 + (5C 2 ioto 7aiCi + <5a 2 C 2 i™, 701 + 5a 2 
a — /3C\C 2 a — /3aia 2 CiC 2 a + /3aia 2 

where a + ftaia 2 =/= 0. 



Remark 2.2. In view of Theorem |2.1| we can obtain by recurrence the distribution of the current angle 
after n steps for the angular random walk 

n 

S n = y arctan Cj 
3=1 

described in the introduction. We have 

n 

S n = /, arctan Cj( aj ..o) = arctan C(a n _ ll o) + arctan C n (a n , 0) 
i=i 

where arctan C(d„_i, 0) is the random variable SVi-i- In particular, if a,j — 1 for j — 1, . . . n, we have the 
following property of the standard Cauchy random variables Cj 

n 

y arctan C j = arctan C 
i=i 

Remark 2.3. A simple byproduct of Theorem |2.1| is that 



„ fOjj 1 f jf, x + y a 1 a 2 dxdy 1 /"2 f* a °1 ta„8 1+ a 2 ta n 8 2 

7T 2 7 K 2 (af + x 2 )(a 2 +y 2 ) ir 2 J_*J_* 

t&M. (2.4) 



e 1+ " 



In (2.4 1 we have used the transformations x = a\ tan#i and y = a 2 tan# 2 . In the special case a\ = a 2 = 1 



the relationship (2.4) yields 



-1*1 = -i r r e ^ rz%Tx dei de 2 = ^ r r e'Wi+^dMfr 

2 /" r 

xcos(/?tan:r)dx. (2.5) 



In the last step of (2.5 1 we have used the transformations #1 + 6> 2 = x and # 2 = y. The integral (2.5 1 
shows that, if (&%, 9 2 ) is uniform in the square 5 = {(#i, 62) ■ — § < \@i\ < § , i = 1, 2}, then the random 
variable W = tan(0i + 2 ) has characteristic function e - '^' because 6i + 2 is uniform and therefore 
W is Cauchy distributed. 

Remark 2.4. It is well-known that for a planar Brownian motion {(B 1 (t),B 2 (t)),t > 0} starting from 
(x,y) the random variable Bi(T y ) possesses Cauchy distribution with parameters (x,y) where 

T y =inf{t>0: B 2 {t) = 0}. 

If the starting points of two planar Brownian motions (B\(t) , B 2 (t)) , for i = 1,2, are located on the y 
axis as in the Figure [4] then we have that 

6 = 6i + 9 2 = arctan B{ (T ai ) + arctan BJ{T a2 ) 



law 



arctan C\ + arctan C 2 
arctan 



1 + ai<2 2 

where C\ and C2 are two independent Cauchy random variables with scale parameters a\ and a 2 respec- 
tively and position parameter equal to zero. Therefore if the starting point of a third Brownian motion 



has coordinates I 0, "i;^ 2 ) then B T a 1 +a 2 I represents its hitting position on the a;-axis. This point 
forms with (0, 1) and the origin a right triangle with an angle equal to 9 = 0i + 82. 
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Figure 4: The hitting position on the x-axis of a planar Brownian motion is Cauchy distributed. In the figure 
the random angles Gi, G2 and G = Gi + G2 are shown. For a = ^+^a ^he right-hand figure shows the hitting 

position B\{T a ). 



3 Non- Centered Cauchy random variables 

For independent Cauchy random variables C\ and C 2l with location parameters 61 and b 2 and scale 



parameters a\ and 02, the random variable U in (2.1l is still Cauchy distributed with both parameters 



affected by the values of the location parameters 61, 62 and the scale parameters oi, 02- 

Theorem 3.1. IfCi, i — 1,2, are two independent, Cauchy random variables with location parameters 
bi and scale parameters a,, then the random variable U is still Cauchy distributed with scale parameter 



ajj 



and position parameter 



bu 



(ai + o 2 )(l + a x a 2 - hb 2 ) + (h + b 2 )(aib 2 + a 2 h) 
(1 + aia 2 - bib 2 ) 2 + (ciib 2 + a 2 bi) 2 



(ai + a 2 )(aib 2 + a 2 bi) - (61 + b 2 )(l + a\a 2 - 6162) 
(1 + aia 2 — b\b 2 ) 2 + (a\b 2 + a 2 b\) 2 



Proof 
We obtain the density function of the random variable U by observing that 



Pr{[7 G dw} ■■ -- E <j Pr <j ^ "t^ E dw 
1 — O1G2 



C, 



and remarking that 
Pr 



C1 + C2 
1 - C 1 C 2 



€ dw 



C 2 = y 



01 



1 + y 2 



l-\-wy *■ 



(l + wy) : 



:dw. 



Therefore 



Pr{[/ G dw} = 



aiQ2 dw 



aia2 dw 

7T 2 

ai<i2 dw 

~^~ 2 



1 + IT 




2 2 (i + ^) 2 (y-& 2 ) 2 + a 2(iy 



1 + y 2 



y - 6i(l + wy)] 2 + a?(l + wy) 2 (2/ - M 2 + ■ 



Ay 



l + y 2 



[y 2 E - 2yF + G] [y 2 - 2yH + K] 



dy, 



where 



E = l + 2b lW + {a{ + b{)w 2 , F = -b x - {a{ + 6f - l)w + b x w 



G = a 2 1 +b 2 1 -2b 1 w + w 2 



H = b 2 



K = a 2 2 + Y 2 . 



We rewrite the integral in (3.1 1 in the following form 

a\CL2 dw 



Pr{[7 G dw} = 



Ay + B 



Cy + D 



y 2 E - 2yF + G y 2 - 2yH + K 



dy 



with 



A + CE = 0, 

-2AH + B- 2CF + DE = 1, 
AK - 2HB + CG- 2DF = 0, 
BK + DG= 1. 



The integrals in (3.2 1 can be worked out by means of the change of variables 



/— F F 2 

yv E = = x\ G . 

Ve V E 



H = x^/K~ H 2 . 



The first integral becomes 



Ay + B 



r°° 1 BE + AF + xA^GE-F 2 
y V E-^ : ) 2 + G-^ V ~ J-°oE (* 2 + 1)VGE^W~ (i 



/El • i' 

and the second one takes the form 

Cy + D 



dy 



D + CH + xCVK - H 2 



(x 2 + 1)VK-H 2 



l-oo (y - B) +K-H 2 
A substantial simplification can be obtained because 

\/GE -F 2 = ai(l + u; 2 ), \J H - K 2 = a 2 . 



dx. 



If we turn back to the distribution (3.2 1, in view of the above calculations, we have that 

BE + AF + xAai(l + w 2 ) 



Pr{[/ G dw} 



a 2 dw 
tt 2 (1 + w 2 )E J_ 



x 2 + 1 



-dx 



+ - 



a x dw f°° D + CH + xCai 



1 + x 2 



dx. 



We observe that, in view of the first equation of (3.3 1, we have that 

aiQ2 dw (A 



E 



C 



1 + x 2 



dx = 



(3.1) 



(3.2) 



(3.3) 



10 



and 



Pt{U € dw} 



BE + AF D + CH \ 
«2 „ n , jv +01 S>dui 

B-CF D + CH'. . 
«2^^— ^v + Oi > du>. 

7T(1 + WT J 7T 



The values of B, C and D can be derived by solving the system (3.3 1, simplified as 



B + 2C(HE -F) + DE = 1, 
-2BH + C{G - EK) - 2DF = 0, 
BK + DG= 1. 

By means of cumbersome computations we arrive at the final result 

Pr{[7 E dw} 

1 (1 + a x a 2 - bib 2 )[w(aib 2 + a 2 &i) + a\ + a 2 ] - (ai&2 + a2^i)[w(l + a x a 2 - b\b 2 ) - {b\ + b 2 ) 
■n [w(l + a x a 2 - b x b 2 ) - (b x + b 2 )] 2 + [w(aib 2 + a 2 6i) + a\ + a 2 ] 2 

(ai+a2)(l+aia2—bib2) + (b 1 +b2)(aib2+a2bi) 
_ (l+a 1 a2 — b 1 b 2 ) 2 + (a 1 b2+a2b 1 ) 2 



(ai+a,2)(aib2+a2bi) — (bi+b2)(l+aia2 — bib2) 
(l+aia 2 — b 1 b2) 2 +{aib2+a 2 b 1 ) 2 



{a 1 +a2){l+aia2—b 1 b2) + {b 1 +b2){aib2+a2bi) 
(l+a 1 a 2 — b 1 b 2 ) 2 + (a 1 b 2 +a 2 b 1 ) 2 



2 • 



Remark 3.1. In view of Theorem |3.1| it is possible to obtain the following particular cases. 
• For a\ = a 2 = 1 and b\ = b 2 = b, we have that 



au 



26 2 + 4 

6 4 + 4 ' 



bu = 



2b :i 



6 4 +4' 



This shows that U has center of symmetry on the positive half-line if b > and on the negative half-line 
if b < 0, therefore the non linear transformation U preserves the sign of the mode. 



• For a\ = a 2 = a and b\ = b 2 = b we have that 

2a(l + a 2 + b 2 ) 
(l + a 2 -6 2 ) 2 + (2a&) 2 ' 



au 



2b(a 2 + b 2 -l) 
(1 + a 2 - b 2 ) 2 + {2abf 



We note that ajj and by depend simultaneously from the scale and location parameters of the random 
variables involved in U. 

4 Continued Fractions 

The property that the reciprocal of a Cauchy random variable has still a Cauchy distribution has a 
number of possible extensions which we deal with in this section. 
We start by considering the sequence 



1 



Vi = 



and show the following theorem. 



V 2 = 



1 



1 + T^ 



v n = 



1 



(4.1) 



i+. 
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Theorem 4.1. The random variables defined in (4-D have Cauchy distribution V n ~ C(a„ b„) where the 
scale parameters a n and the location parameters b n satisfy the recursive relationships 



Qn+l — 



(i + M 2 + < 



b n +l 



(i + M 2 + «r 



ra = l,2,... 



n = l,2,.. 



(4.2) 
(4.3) 



Proof 
Let us assume that V n possesses Cauchy density with parameters a n and b n , therefore V n +i writes 



Pr{V„ +1 <„} = Pr{^<„}=P r { TT ^ 



After some computations the density of V n +i can be written as 
fv„+i 0) = ■ 



(l+6„) 2 +< 



(l + b „)2 +a 2 



< V 



v e 



[(l+b„) 2 +a^] 2 



It can be directly ascertained that V\ possesses Cauchy distribution with parameters a\ = 1/2 and 
6i = 1/2. ■ 



Remark 4.1. We have evaluated the following table of parameters a n and b n : 

n 1 2 3 ... 10 2 






1/2 
1/2 



1/5 
3/5 



1/13 
8/13 



5.77e- 4 ^ 
0.618034 



For n — 1,2,3 we can observe that the scale parameters a n coincide with the inverse of the odd-indexed 
Fibonacci numbers while the sequence b n has the numerators coinciding with the even-indexed Fibonacci 
numbers and the denominators correspond to the odd-indexed Fibonacci numbers. 
In light of the previous considerations we can show that for n > 1 



b n = 



F2n 

F 



1 



2n+l 



F, 



2n+l 



(4.4) 



where F n . n > is the Fibonacci sequence. Recalling that the Fibonacci numbers admit the following 
representation (it can be easily checked by induction) 



b n - (1 



75 



(4.5) 



where <\> = ^ is the golden ratio, we now prove that if a n and b n have the representation (4.4 1, then 



also a n+ i and b n+ i can be expressed in the same form. From (4.2 1 and (4.3 1 we have 

F 2 „ 



F-2„ 



1 



->n+l 



F 2n 
F 2 „+i 

F2n+2 



1 



-F2ra+2^2n+l 

~F2 i 1 

r 2n+2 ' 1 



k2n+l 



-(1- 



\2n+l 



F2n 



^2 



04n+4 + (1 _ 0)4«+4 _ 2(f) 2n + 2 (1 - (j)) 2n + 2 + 5 

^2„+l _(!_ 0)2«+l 
[02™ + ! - (1 - 0)2»+ 1 ] [0 2 »+3 - (1 - 0)2n+3] 



>/5 

75 
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F 2n + 2 
F 2n +3 



Similar calculations prove that a„ +1 = jA — . In view of representation (4.4) and (4.5), it is easy to show 



that 



lim b n — lim = lim 

n— foo n— >oc x*2n-l-l n 



m 



*-(i-«( i ?) 



2» 



Otherwise, observing that the sequence b n , n > 1 is increasing, because 
bn+i = F 2n+2 F 2n+1 = 4 "+ 3 + (1 - 0) 4 »+ 3 + 1 

b„ F 2n+3 F 2n ^+3 + (1 _ 0)4n+3 _ 4 



lim a n = lim 

n— »-oo n— foo .FOn-l-l 



> 1 



and taking the limits in (4.2 1 and (4.3 1 we have that 



L = 



L 



{1+H) 2 + L 2 ' 



iJ = 



H+l 



{1 + H) 2 + L 2 ' 



(4.6) 



where iJ = limn^oo a n and i = linin^oo b n . From the relationships in (4.6) we derive the equality 

L L 



H H + l 
that implies L = 0. In fact, for L ^ 0, we arrive at the absurd that H = H + 1. Substituting L = in 

if X 



the second formula of (4.6 1 we obtain 



H + l' 

since H satisfies the algebraic equation H 2 + H — 1 = it follows that H — </> — 1 where is the golden 
ratio (see Figure 5). 



Remark 4.2. A slightly more general case concerns the sequence 
1 1 ... 1 



Wi 



ci + rfiC( a0i6o ) ci + a a di + b Q di 



W 2 



c 2 + d 2 Wi 



W 3 



1 



c 3 + d 3 W 2 



By performing calculations similar to those of Theorem |4.1| we have that W\ has Cauchy distribution 
with scale parameter a\ and position parameter b\ such that 



«i 



d\a 



(ci +rfi6 ) 2 + rf?c 



61 = 



a +d 1 b 



( Cl + dibo) 2 + d 2 a 2 ' 



Similarly, if W n ~ C(a n ,&„)> than W n+ i - C'(a„ +1 ,fc„ +1 ) where 

«n+l a n 



flra+1 



(cro+i + d ra +i&„) 2 + <^ +1 a 2 ' 



y n+l 



c n +i + d n+ ib n 



(c„+i + d„ + i6„) 2 + d 2 +1 a 2 



(4.7) 



for every n > 2. The sequences in (4.7 1 for c n = d n = 1 coincide with (4.2) and (4.3 1. 

Another sequence of continued fractions involving the Cauchy distribution is the following one 

1 1 „ 1 



Ui 



1 + C 2 



U 2 



1 



u„ 



1+c 2 



1 



(4.8) 



^■■■i+W 



It is well-known that the random variable U± possesses the arcsin law. Unlike the sequence V n studied 
above the sequence U n , n > 1, has a density structure changing with n. Some calculations are sufficient 
to show that U\, U 2l U 3l U4 have density, respectively equal to 



A/i ( u ) 



1 



w-\/u(l — u) 



< u < 1, 
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Figure 5: In the left figure the densities of the Cauchy random variables Vi, V2, V3 and V4 are shown. In the 
right figure the densities of U\, U2, U3 and U4 are plotted. 



fu 2 («) 
fu 3 {u) 

fu 4 (u) 



Tiuyf {I - u){2u - 1)' 
1 



tt(1 - u)>/(2u-l)(2-3w)' 
1 



2 <u<l, 

1 2 

2 <U< 3' 
3 2 

5 3 



tt(2u - l) v /(2-3u)(5u-3) : 
The general result concerning U n is stated in the next theorem. 
Theorem 4.2. For every n > 1 i/ie distribution of the random variable U n is given by 

Pr{U n e du} 



wht 



7r[(-l)»+i a „ + (-1)"/M ^(-l)»/3 n + (-l)™+i(a n + /3 n ) u 

x — — du, 

V(-1)" +1 K + /?„) + (-l) n K + 2/3„)u 

{ _ 1)n ^n+j3n < < ( _ 1} n_^_ 

a n + 2/} n a„ + p n 



and a n ,/3 n € N satisfy the recursive relationships a n — f$ n -x, j3„ — a n -i + /?«— i- 
Proof 



From (4.8 1 we have that 



fn+l — 



1 



1 + J7„ 



then proceeding by induction, i.e. assuming that U n has distribution (4.9), we obtain that 

d 



Pr{(7 n+1 G du} 



1-ul d r h ( a «'M 

. -Pr <j U n > }du=— / Pr{U n E du} 

du u du /1-1 



1 1 



1 



1 



(4.9) 



vr U * [(-Ijn+lo,, + (-!)»&(*=*)] ^ + (_i)n + i (an + ^)(i^) 
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I 

i 1 



^(-l)«+ 1 (a„ + p n ) + (-!)«(«„ + 2/3„)(i^) 



tt[(-1)»/3„ + (-l)»+ 1 (a n + /3„)u] ^(-lJn+^Qn + /3„) + (-l)«(a n + 2/3> 

x — du. 

v/(-l)"K + 2/3„) + (-l)»+i(2a n + 3/3 n )u 

In the first integral the function h(a n ,j3 n ) represents the right boundary of the support of U n . We con- 
clude that t/,1+1 possesses distribution fi~9| by taking a n+ \ = /3 n and (3 n +i = cn n + f3 n . ■ 



Remark 4.3. The sequence /3„ is a Fibonacci sequence since we have that /3 n = f3 n ~i + On-i = Pn-i 
Ai-2- We note that the sequence of coefficients a n and j3 n are such that 

1- a n+l y Pn , 1 

lim — = lim — = <p — 1 . 

n-J-oo ^jj+j n^oo fj n+1 



On the base of arguments similar to those of Remark |4.1| it is possible to show that the sequence U n , 
n > 1, converges in distribution to <f> — 1. In this case the upper and lower bounds of the domain of 
definition of the densities fu n (u), n > 1 are expressed as ratios of Fibonacci numbers (see Figure 5). 
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